Nucleon-nucleon charge symmetry breaking and the dd — > a7r° reaction 
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We show that using parameters consistent with the charge symmetry violating difference between 
the strong nn and pp scattering lengths provides significant constraints on the amplitude for the 



dd 



air reaction. 
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The concepts of charge independence and charge sym- 
metry provide powerful tools in organizing the multi- 
plet structure of hadrons and nuclei. These symme- 
tries are not perfect; diverse small but interesting vio- 
lations have been discovered [l|, Q • Our concern here 
is with the breaking of charge symmetry (CS). This 
symmetry is defined as invariance under a rotation by 
180° around the 2-axis in isospin space. In quantum 
chromodynamics (QCD), CS implies that dynamics 
are invariant under the exchange of the up and down 
quarks However, since the up and down quarks 
do have different masses (m u £ m d ) [3], the QCD 
Lagrangian is not charge symmetric. This symmetry 
violation is called charge symmetry breaking (CSB). 
The different electromagnetic interactions of the up 
and down quarks also cause CSB as well as the break- 
ing of charge independence. Observing the effects of 
CSB interactions therefore provides a probe ofm u and 
nid, once the electromagnetic interactions are treated. 

It has long been known that CSB is violated in the 1 S'o 
state of nucleon-nucleon scattering, with a pp — a nn = 
Aa = 1.5 ± 0.5 fm [l|, where a denotes the scatter- 
ing length. The nn interaction is more attractive than 
the pp. There are a variety of explanations for this us- 
ing meson exchange mechanisms [1, [f| . Each of these 
mechanisms involving the strong interaction can be 
traced to the mass difference between the up and down 
quarks. Nucleon-nucleon potentials that are consis- 
tent with this scattering length difference are success- 
ful in reproducing (along with electromagnetic effects) 
the measured binding energy differences between mir- 
ror nuclei [5j, |6| . It is interesting to search for further 
manifestations of the up-down quark mass difference. 

Two exciting observations of CSB in experiments in- 
volving the production of neutral pions have stirred 
interest in this subject. CSB was observed in the re- 
action np -> cZtt at TRIUMF by measuring the CSB 
forward-backward asymmetry of the differential cross 
section as A {h = [17.2 ± 8(stat) ± 5.5(sys)] x KT 4 



Q- Furthermore, the final experiment at the IUCF 
Cooler ring reported a relatively large dd — > air cross 
section (a = 12.7 ± 2.2 pb at T d = 228.5 MeV and 
15.1 ± 3.1 pb at 231.8 MeV) %. The dd -> an reac- 
tion violates CS since the deuterons and the a-particle 
are self-conjugate under the CS operator, with a pos- 
itive eigenvalue, while the neutral pion wave function 
changes sign. 

The study of CSB n° production reactions presents an 
exciting new opportunity to determine the influence 
of quark masses in nuclear physics, and to use effec- 
tive field theory (EFT) to improve understanding of 
how QCD works This is because chiral symmetry 
of QCD determines the form of pionic interactions. 
Electromagnetic CSB is typically of the same order of 
magnitude as the strong one, and also can be handled 
using EFT. The EFT for the Standard Model at mo- 
menta comparable to the pion mass, Q ~ m T , is chiral 
perturbation theory (xPT) M. This EFT has been ex- 
tended to pion production [TfljlDJ QJ, GJ, 03, 

where 

typical momenta are Q ~ y/m^M, with M the nu- 
cleon mass. EFT using the operators of [15] was used 
to correctly predict the sign of the forward-backward 
asymmetry in np — > dw [16| . 

The purpose of the present note is to use informa- 
tion regarding CSB in the nucleon-nucleon system to 
constrain or inform the calculations of the dd — > an 
reaction. We begin by describing the pion production 
calculations and then show how the nucleon-nucleon 
CSB is relevant for this calculation. We constrain 
the parameters by the requirement that NN CSB is 
consistent with observation. It turns out that this 
constraint reduces considerably the uncertainty in the 
predictions for the dd — > air reaction. 

Attempts to understand the dd — ► an reaction be- 
gan with a survey of various mechanism using ini- 
tial state plane wave functions and simplified final 
state wave functions [ll| . Next, recent significant ad- 
vances in four-body theory [l7], EH were used to in- 
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FIG. 1: Diagrams of np — » d7r°; the solid circle indicates 
CSB. 

elude the effects of deuteron-deuteron interactions in 
the initial state, and to use bound-state wave func- 
tions with realistic two- and three- nuclcon interactions 
The resulting calculations are hybrid: the CSB 
pion-production operators are constructed using EFT, 
but the strong interactions are not. The result was 
that a cross section of the experimentally measured 
size could be obtained using leading order (LO) and 
next-to-next-to LO (NNLO) pion production opera- 
tors. However, not all NNLO diagrams were included. 
A complete analysis would require a careful treatment 
of loop diagrams. 

The calculation made use of a variety of CSB mech- 
anism which we will briefly review now. At formally 
leading order, there is only one contribution, repre- 
sented by Fig.[Th.: pion rescattering in which the CSB 
occurs through the seagull pion-nucleon terms linked 
to the nucleon-mass splitting. 

There is no next-to-leading order (NLO) contribution. 
At NNLO, suppressed by 0(m v /M), there exists a 
recoil correction of the LO term, labelled by Ai lcc . 
This is determined by the same parameters as the 
LO term. At the same order there are other opera- 
tors. A term arises in which a one-body CSB operator 
(oc /?i + fls) is sandwiched between initial and final 
state wave functions, as illustrated in, e.g., Fig. [TJd. 
We refer to this as the one-body term (A^ib). The 
terms (3\ = 0(em^/M 2 ) and /?3 = 0(a/ir) arise from, 
respectively, the quark-mass-difference and electro- 
magnetic contributions to the isospin-violating pion- 
nucleon coupling. Neither (3\ nor /?3 can be extracted 
from experiment yet. These terms were estimated by 
modeling [20j f3\ by ir-rj mixing, see Fig. \Tjp, 

(3 1 =g 7 ,(ir°\H\ri)/m% (1) 

where (iv \H\rj) — —4200 MeV 2 is the 7r-?7-mixing ma- 
trix element [2l| . and g v = griNNf^/M the 77-nucleon 
coupling constant. Nuclcon- nuclcon clastic scattering 
data show little sensitivity to n exchange and high- 
accuracy fits can be achieved [22| using g^ NN /^ = 0. 
Indeed, the possibility of a vanishing coupling con- 
stant had been raised earlier by the detailed analysis 
of NN total cross sections and pp data using disper- 
sion relations [23[ resulting in g\ NN j^ ~ 0. Also, in 



the Bonn full model [24j it is found that g^ NN /4:TT w 
is consistent with the NN scattering data. A value of 
dr/NN/^ 77 = 0-51 was use d m Ref- 03 ■ 

The effects of electromagnetic interactions as well as 
strong CSB were included in computing the a particle 
wave functions. These interactions generate a small 
isospin T = 1 component of the wave function that en- 
ables a non-zero contribution of charge-symmetry con- 
serving (CSC) production operators. The one-body 
operator was used to generate n° production . 
A number of other CSB mechanisms enter at N 3 LO or 
higher, including loop diagrams, and short-range in- 
teractions. The lowest order, where four-nucleon con- 
tact interactions start to contribute, is N 4 LO. To es- 
timate their strength, Ref. [TTf| evaluated certain tree- 
level contributions as indicated by Figs. QJ;. This fig- 
ure represents the exchange of heavy mesons (a, u>, p) 
via a Z-graph mechanism, with tt-t] mixing to generate 
CSB at pion emission (Af CT , Af w , and M. p ). Another 
Z-graph (labeled as A4 puj ) arises in which the CSB 
occurs in the heavy-meson exchange via p-u mixing 
along with strong pion emission at the vertex. 
The Z-graphs are believed to be important, because 
their inclusion lead to a quantitative description of the 
total cross section for the reaction pp — > ppir near 
threshold [25| . The results fl9ll u se the coupling con- 
stants and parameters of Ref. [IlJ, see their Table I. It 
was found that the Z-graphs give unexpectedly large 
contributions, especially the p-uj exchange operator 
that add constructively and overwhelm the one-body 
term. This model of resonance saturation, gives re- 
sults in vast disagreement with the power counting 
and therefore needs reassessment. Here we re- asses 
the coupling constants used by [l9[. This is only a 
first step, because it is also necessary to justify the 
use of Z-graphs in resonance saturation, which would 
require further calculations. 

The various mechanisms generate pion-production 
kernels that are sandwiched between final and initial 
state wave functions to provide a transition matrix el- 
ement M. for Td— 228.5 MeV. These matrix elements 
are given in Table 1 of [19j . The transition amplitude 
can be written as 

M = MpE + M.cc + Mib + Ma+Mp + M,, 

+M pu1 +Mwf, (2) 

where the pion exchange term A^pe, its recoil cor- 
rection A^rec, and the effects of CSB in the a wave 
function A4\vf are independent of g v and (3\. The 
one-body term and sigma and rho exchange terms in- 
volving the Z-graphs, A'lib + M. a + M. p , are propor- 
tional to pi. The terms M.^ + A4 pu >, which arise from 
omega and rho-omega exchange, are proportional to 
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(3\ and g^, the strong w-nucleon coupling constant. 
Given the contributions to M. expressed in Tabic 1 of 
in units of 10~ 4 fm~ 2 , the cross section can be 
written as 



a = 4.303 pb \M x 10 4 fm 2 



(3) 



Results were obtained [19| using either the Argonne 
Vis (AV18) [H] or CD-Bonn (CDB) two-nucleon 
potentials combined with a properly adjusted Tucson- 
Melbourne (TM99) [13] three-nucleon force. The com- 
bination guarantees that the a particle binding energy 
is reproduced with high accuracy. 
We now turn to CSB in the nucleon-nucleon (NN) sys- 
tem due to meson-mixing. Our CSB NN calculation 
is described in Ref. [5[. Since we have to use meson- 
nucleon form factors in our NN scattering calculations, 
while in Refs. [Til . Il9j no such form factors are ap- 
plied, we explain the definition of meson-nucleon cou- 
pling constants in conjunction with form factors. For 
this, we define a coupling constant as a function of the 
momentum-transfer t by 



g a (i) = F a (t)g a (t = to 2 ) 



with 



F a {t) 



A 2 



(4) 



(5) 



where a stands for any meson, m a denotes the meson 
mass and A Q the so-called cutoff mass. We use A PiW = 
1400 MeV. In the relativistic three-dimensional the- 
ory used for CD-Bonn, the momentum transferred be- 
tween the two nucleons is t — — (q' — q) 2 with q and q' 
the center-of-mass nucleon three-momenta before and 
after scattering. Coupling constants used in theories 
without form factors should be compared to coupling 
constants at t = 0, i.e., g a (0), of a theory with form 
factors. Thus, in our CSB NN calculations, we use 
the coupling constants given in Table I of Ref. 
and identify them with <? Q (0), except for the omega- 
nucleon coupling where we use 



g£(0) 

47T 



5.0 



(G) 



instead of the 10.6 applied in Ref. [Tl|. Our lower 
value for the omega coupling is more consistent with 
a fixed- s dispersion relations analysis by Hamilton and 
Oades [28| in which a value of 5.7 ± 2.0 was obtained. 
Moreover, SU(3)a avor implies g u = 3g p which, for 
the rho-coupling used, yields g 2 /47r — 3.9. Apply- 
ing the omega coupling constant stated in Eq. (|6|) and 
the other meson parameters as given in Table I of 



Ref. [1111 and using (p°\H\: 



-4300 MeV 2 for the 



p-u-mixing matrix element 21|, we find a CSB con- 
tribution to the 1 So scattering length difference from 
p — uj mixing of 



1.45 fm. 



(7) 



Obviously, this term entirely accounts for the observed 
CSB scattering length difference. We have also calcu- 
lated T] — 7T° mixing to find 



Aet„ 



0.33 fm 



(8) 



using the same strong coupling constants and mix- 
ing matrix element as in jl9j |. Since the puj mixing 
explains all CSB, this might be considered as yet an- 
other argument that g v = 0. 

To explore a variety of possibilities, we re- write the 
results of 113 as 



M 



CDB 



-2.75 + 3. li + (-3.4 + 2. 82i) 



g 2 /4^ 
.51 



+(-.53 + .44i) 
+(-1.32 + 1.51i 



.51 10.6 



10.6 



.51 - .13z) 10 -4 fm" 



A4avis = -1.65 + 1.58?+ (-2.44 + 2.2H) 



(9) 



g 2 /4^ 
.51 



-(-.35 + .34i) 



+(-0.84+1.07i) 



g'/47r g 2 /47r 
.51 10.6 



g 2 /4?T 



10.6 



.41 - .14i) 10" 4 fm" 



(10) 



depending on the potential used. Ref. [l9( observed 
that this model dependence visibly influences the cross 
section result, requiring a more consistent treatment 
of the NN interaction and production operator in the 
future. 

As a first step, we note that using the parameters of 
11011 leads to the results 



A4cdb = I 

A4 A V18 = I 

with cross sections 



-7.49 + 7.74i)10" 4 fnr 2 , 
-4.87 + 5.06i)10- 4 fm~ 2 , 



0CDB 
CAV18 



499 pb, 
212 pb. 



(11) 
(12) 



(13) 
(14) 
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Next we replace g 2 / Air by the value stated in Eq. © , 
which yields 

A^cdb = (-6.80 + 7.03i)lCT 4 far 2 , (15) 
Mavis = (-4.42 + 4.54i)l(T 4 far 2 . (16) 

This would lead to about a 20% reduction of the cross 
section. However, using 7r?7 mixing along with pcu mix- 
ing overestimates Aa. Therefore, it is reasonable to 
explore the consequences of using g v = 0, which re- 
sults in 

Mcbb = (-3.15 + 4.00i)10~ 4 far 2 , (17) 
Mavis - (-1-82 + 2.17i)10~ 4 far 2 , (18) 

with cross sections 

acDB - 111.5 pb, (19) 
cTAvis = 34.5 pb. (20) 

The cross section for AV18 is now only a factor of 2 
or so bigger than the data, while CD-Bonn is off by a 
factor of about 7. 



In summary, past theoretical work (l9| on the cross 
section of the dd — ► air reaction at 228.5 MeV was 
plagued by the problem that the predictions were off 
by factors between 15 and 30. In this note, we have 
shown that constraining the coupling constants in- 
volved by the requirement that the CSB in the 1 So 
NN scattering length is correctly reproduced reduces 
the over-prediction to just a factor of about 2 (using 
the AV18 potential). In relative terms, this is substan- 
tial progress in understanding the dd — > air reaction. 
However, significant differences between the use of the 
AV18 and CD-Bonn potentials remains, signaling that 
a deeper understanding is needed. 
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